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ABSTRACT 

The recent discovery of luminous quasars up to a redshift z = 6.43 has renewed interest in the formation of 
black holes massive enough to power quasars. If black holes grow by Eddington-limited gas accretion with a 
radiative efficiency of at least 10%, the time required to grow from a stellar black hole to ~ K^Mq is ^ 10^ 
years, close to the age of the universe at z = 6. Black hole mergers may accelerate the rate of mass growth, but 
can also completely eject black holes from halo centers owing to the gravitational wave recoil effect. Recently, 
Haiman concluded that black hole ejections likely do not allow black holes to grow to ~ lO^M© by z = 6.43. We 
reexamine this problem and show that, by using a different halo escape velocity, accounting for the dependence 
of the recoil velocity on the black hole binary mass ratio and spins, and allowing seed black holes to form in all 
halos down to virial temperatures of 2000 K, black hole masses may reach ^ lO^M© as early as z = 9 starting 
from stellar seeds, without super-Eddington accretion. In this particular case, we find that these massive black 
holes form from the merger of '-^ 10"* stellar black holes formed in low-mass halos at z '^ 20, which must all 
grow close to the maximum Eddington rate over most of the time available from their birth to z ~ 6. The 
alternative is that black holes can grow more rapidly by super-Eddington accretion. 

Subject headings: black hole physics — cosmology: theory — galaxies: nuclei — gravitation — gravitational 
waves 



1. INTRODUCTION 

Highly luminous quasars at high redshift require the for- 
mation of black holes (BH) with a mass M large enough 
to make the observed luminosity be less than the Edding- 
ton limit. This may pose a problem for cold dark matter 
(CDM) models where massiv e halos that can harbor BHs do 
not form until late epochs (.Efstathiou & Reeslll988h . Cur- 
rently, the h ighest redshift qu asar known is SDSS 1 148H-325 1 , 
at z = 6.43 (iFan et alJl2003l) . Its luminosity L, if equated to 
the Eddington value Ledd = 4 x \Q'^Lq{M /Mq), results in a 
minim um BH mass M = 4 x 1O''M0 (Fan et al. 2003; Haiman 
l2004h . The quasar luminosity inferred from the observed flux 
is not likely to be g reatly overestirnated due to gravi tational 
lensing or beaming ( Haiman & Cenl2002tlwrilott et alJ2003t 
Eeeton et al. 2004). 

If BHs grew in mass by accretion through a standard thin 
disk with radiative efficiency e ^^ 0. 1 , then a BH of initial mass 
Mq formed at time fo accreting continuously up to time t with 
luminosity L will reach a final mass M = Mq exp[(f-fo)/fsai], 
where fsai is the Salpeter time (Salpeter 1964), 



radiative efficiency is e < 0.1 (e.g., Ohsuga et al. /OOz), or 
the luminos ity is grea ter than the Eddington luminosity with 
e = 0.1 ( Ru szkowski & Begelman 2003). (c) Finally, several 
BHs may have merged to form a more massive BH. 

In this paper we consider the latter possibility in the con- 
text of the CDM model, where dark matter halos merge hi- 
erarchically. If the earliest BHs arose from the first stars 
forming in the center of halos, then h alo mergers can lead t o 
the merger of their central BHs (e.g., IVolonteri et alJl2003l) . 
after dynamical friction leads to the formation of a binary 
BH which can then lose orbital energy by interacting with 
stars (Begelman, Blandford & Rees 1980) or via gas dissipa- 
tion dGould & R ix 2000), until the emission of gravitational 
waves takes over If many BHs growing from initial stellar 
seeds of mass Mq. made at time fo, continuously accrete un- 
til time t (as they merge with each other at any intermediate 
times), the final BH can reach a mass ( Haiman & Loeb 2001 ), 
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Note that the factor 1 - e in the denominator is to account for 
the fact that only a fraction 1 - e of the accreted mass is added 
to the BH. Since the age of the universe at z = 6 is ~ 10^ years, 
a BH made by one of the first stars formed in the universe with 
Mq ^ IOMq would need to be accreting near the Eddington 
rate for almost all the time available before z = 6 in order to 
increase its mass by a factor ^ 10^ to power the observed 
quasars (corresponding to ^ 20 e-folding times). 

There are several ways to ease this requirement and allow 
BHs to grow to a large mass in less time: (a) BHs may have 
started at a high er initial mass from the collapse of supermas- 
sive stars (e.g., iCarr. Bond & Arnett.1984.) . (b) There may 
be super-Eddington accretion with fsai < 4 x 10^ years if the 



In the last stage of the merger of BHs, the emission of 
gravitational waves (GW) can give a net momentu m to the 
resulting BH large enough to eject it from its halo (TitcheM 
1983; Kidder 1995; Favataetal. 2004; Madau et al. 200^ 
iMerritt et al J 12004ft . Then the halo center is left empty and 
the process of BH growth must start over again. Recently, 
iHaimanl ( 12004ft presented a calculation of the growth of BH 
masses assuming that they start growing from stellar seeds 
only when the halo has a velocity dispersion a > Vkwk/'^, 
where Vkici< is a fixed kick velocity given to the BHs by the 
final GW burst. Haiman (2004) concluded that for ua,ca > 
64 kms~', BHs would not be able to grow to large enough 
masses without super-Eddington accretion. 

We show in this paper that when the adequate kick velocity 
as a function of the mass ratio and spins of the two merging 
BHs is used, and a halo escape velocity based on an isother- 
mal profile (higher than assumed by Haiman) is used, ejec- 
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Fig. 1 . — Mass of the main BH progenitor as a function of redshift ending 
at z = 6.43, according to four different estimates of GW recoil effect (shown 
as the four lower curves). The upper curve (thick solid line) shows the mass 
of the main halo progenitor Seed BHs ai'e started at IQMq, and the maxi- 
mum BH mass allowed is 10"^ of the halo mass. For reference, GW recoil 
velocities are of order 1000, 1000, 400, and 50 kms"' for Fitchett, Kidder, 
Favata upper and lower models, respectively. 



tio ns by GW emis sion result in a smaller reduction than found 
by Haiman ( 2004) of the final black hole mass that is attained. 
As a result, we shall show that BH masses of IQ'^Mq can be 
reached up to z ~ 9 without super-Eddington acc retion, eyen 
for the highest plausible kick velocities found by Favat a et aO 
(|2004). Our model for the BH evolution and growth is de- 
scribed in § 2. The main results are presented in § 3, and we 
discuss the implications in § 4. 

2. BLACK HOLE EVOLUTION MODEL 

We model the evolution of dark matter halos using the 
merger tree of Cole et al. (2000) ba sed on the ext e nded Press- 
Schechter formahsm (.Press & Schechteri 119741 iBond et aTl 
fT99ll iLacev & Colelll993l) We use the flat ACDM power 
spectrum ("Eisenstein & Hu"1999') with matter density r2„, = 
0.3, baryon density fth = 0.04, present Hubble constant Ho = 
70 kms~' Mpc~', power spectrum normalization ag = 0.9, 
and primordial index n = 1, consistent with the WMAP re- 
sults (Bennett et al. 2003; Spergel et al. 2003). We start the 
merger tree at z = 6.43 with a final halo mass Mf, find- 
ing all the progenitor halos at higher redshift. The mass 
Mf is determined by requiring the quasar number density 
hq = 2.7 X IQ-'^ih-^ Mpc)3 ("F an et alJ20()l to be equal to the 
number density of halos of mass Mf at z = 6.43. This gives 
the maximum possible mass of the host halo of the quasar 
consistent with the observed quasar abundance (for example, 
if quasars have a short duty cycle the host halo abundance 
must be higher than the quasar abundance, so the halo mass 
must be lower) . The resulting halo mass, Mf ~ 8.5 x 1O'^M0 
(jHaiman 2004), requires a 5.2a fluctuation at z = 6.43. 

The first BHs are assumed to form in halos that have 
reached a virial temperature T = 2000 K, when molecular hy- 



drogen cooling indu ces the formation of the first stars (e.g., 
1 120031) . A constant initial BH mass Mq is as- 

The 



Yoshida et all 

sumed in every halo reaching this virial temperature, 
mass corresponding to this virial temperature, Mres(z), is cho- 
sen as the resolution of the merger tree. In order to accurately 
follow the appearance and growth of halos down to mass M^.,, 
we compute accretion into halos of mass M/, from smaller ha- 
los down to a minimum mass 10~''M/, or Mres, whichever is 
smallest. We compute the virial temperature of halos, related 
to velocity dispersion as T = fia^/k (where k is Boltzmann's 
constant and /i is the mean particle mass), using the fitting 
formula of Brvan & Norman ( 1998) for the mean overdensity 
of virialized halos A^, and the equations a^ = GM/,/(27?/,), 
Ml, = {4tt /3)RIpc{z)Ac(z), where R/, and pc are the radius of 
virialization of the halo and the critical density at a given red- 
shift z, respectively. 

We assume that all the BHs grow continuously by gas ac- 
cretion on the timescale fsai obtained with e/(l -e) = 0.1 and 
L = LEdd, after their host halo has reached a minimum veloc- 
ity dispersion cr,„,„. We shall show some results for ct,,,,-,, = 0, 
in which case BHs start growing as soon as they are formed. 
Alternatively, the effects of slow gas cooling, stellar winds, 
or supernovae may prevent any gas concentrating in the halo 
center and being captured by the BH for accretion when 
<T < amin- The BH growth by gas accretion continues until 



the BH mass reaches a fraction 10 of its halo mass, at which 
point we simply set M = IQT^Mk- 

The last step is to check if the BH merger results in an ejec- 
tion from the host halo. We consider four estimates of the 
GW recoil velocity. The first is a quasi-Newtonian calcula- 
tion of non-spinning BHs by Fitchett ( 1983). Second. Kidden 
( 1995) added a post-Newtonian spin-orbit correcti on to Fitch- 
ett's work that depends on the BH spins. Lastly, Fava ta et all 
(2004) and Merritt et al. (2004) obtained a new estimate using 
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Fig. 2. — Composition of the final BH. Upper panel: Number of seed 
BHs that merged into the final BH. as a function of their formation redshift. 
L{)wer panel: Contribution to the final mass of the BH from the BH seeds as 
a function of their formation redshift. 



BH perturbation theory, and they provided a lower and upper 
limit on the GW kick velocity comprising the uncertainty due 
to the final plunging state which gives the dominant contri- 
bution (see their eqs. [1], [2]). We use these limits as two 
separate estimates. 

The mass acccretion histories of BHs are calculated fol- 
lowing the merger trees in time and computing the GW re- 
coil velocity every time two halos merge (and by assumption, 
their central BHs merge). The dimensionless spin parame- 
ters of merging BHs are uniformly chosen from antiparal- 
lel (fl = -1) to parallel (a = 1) in relation to the orbital an- 
gular momentum, to calculate the effective spin parameters 
iDamouiil200U iDamour & Buonannoll200(il) . We then com- 
pare the GW kick velocity with the escape velocity of the 
halo to decide if the merged BH is ejected. We compute the 
escape velocity assuming the halo has an isothermal profile, 
with density p oc r~^ down to the radius of the BH zone of 
influence: R~ = GM /2a^. Therefore, the escape speed at R. 
is Vesc = 2CT-\/l+ln(7?/,/7?;), or about 5 a when M = \Q~^Mh. 
Note that Haiman (2004) assumed instead that Vesc = 2a, and 
that no BHs start growing until the halo escape velocity has 
exceeded an assumed, fixed GW kick velocity. The exact 
value of the escape velocity depends on the halo density pro- 
file, but it would not be much lower than we assume for a pro- 
file similar to observed galaxies. Gravitational waves emis- 
sion also implies that a fraction of the mass in BHs is lost as 
the energy of the waves. T his loss of energy does no t greatly 
reduce the final BH mass (iMenou & Haimanll2004h . and we 
neglect it in this paper. 

3. RESULTS 

Figure [2 shows the mass history of the BH in the z = 6.43 
quasar SDSS 1 148H-525 1 in our models. The upper thick solid 
line shows the mass of the main halo progenitor in one real- 
ization of the merger tree, obtained by always choosing the 
branch of the most massive progenitor at every merger The 
lower four lines show the mass of the BH in this main halo 
progenitor, according to our four different prescriptions for 
the GW kick velocity. The mass of the seed BHs is fixed at 
Mo = lOM©, and they start growing by gas accretion immedi- 
ately as they are formed, over a time fsai = 4 x 10^ years. At 
every halo merger, the mass of the two BHs is added when 
the GW kick velocity is smaller than the escape velocity. This 
results in the sudden BH mass increases seen in the figure. 
Growth by gas accretion is assumed to continue immediately 
after mergers. When the GW kick velocity exceeds the escape 
velocity, the BH is removed and replaced immediately with a 
new seed with Mq = lOM©. The figure shows that BHs are of- 
ten kicked out at high redshift when they reside in low-mass 
halos, but this stops when the halo reaches a high enough es- 
cape velocity. Large, sudden mass increases (for example, at 
z ~ 13 for the Kidder curve) are due to additions of BHs that 
have grown to high mass in a merging subhalo. Note that 
the Favata lower limit curve reaches a mass 10"^ M/, (at which 
point the BH is not allowed to grow further) at z ~ 10. 

Figure|2]shows a histogram of the number of BHs that have 
merged into the final BH at z = 6.43, as a function of their for- 
mation redshift (top panel), and their contribution to the total 
mass of the final BH [bottom panel). In other words, the num- 
ber shown in the top panel is the number of BHs that we need 
to add over in equation Q from every interval of BH forma- 
tion redshift, and the bottom panel shows their contribution to 
the final mass including the factor gf'"'", V'sij xhe figure shows 
that the final mass attained by the BH at z = 6.43 is the result 




9 10 15 20 30 40 

redshift (1+z) 



9 10 15 20 30 40 

redshift (1+z) 



Fig. 3. — BH mass histories as in Fig. 1, clianging model pai'ameters for 
mergers and growth by gas accretion. Upper left panel: the seed BH mass is 
increased from WMq to lOOA/0. Lower left panel: Growth by gas accretion 
stalls only once the host halo has reached a minimum velocity dispersion 
o"nrai = 15 kms"' . Upper right panel: BH mergers occur only when merging 
halos have a mass ratio higher than q = 0.\. Lower right panel: The minimum 
velocity dispersion is increased to cr,,,,,, = 50 kms"'. 



of the merger of many small BHs: the formation redshift bins 
with the highest contribution to the final BH mass typically 
contain lO'* to 10^ seed black holes that merged. Hence, the 
~ 10** mass growth factor from the initial seed mass of lOM© 
to the final BH mass of ~ lO'^M© originates in about equal 
parts from the growth factor e('"'»,)/'sai of each individual BH, 
and from the number of similar terms that are added together 
in equation (|2j from all the BHs that have merged. 

In summary, these results show that in the optimistic model 
we have assumed for BH growth, it is possible to reach a 
mass as high as that inferred for the BH in the q uasar SPSS 
1148H-5251 at a redshift of up to z ~ 10 for the Favata etal J 
( 2004) lower limit of the GW kick velocity. Even for the high- 
est GW kick velocities the same BH mass can still be reached 
atz~8. 

In Figure |3] we examine a number of variations on our ba- 
sic model. Increasing the seed BH mass to lOOM© {upper 
left panel) causes a small increase of the redshift at which a 
mass of ^ IO^Mq can be achieved. The models in the lower 
left panel show what happens when growth by gas accretion 
is allowed only in halos with cr„„„ > 15 kms~', since it may 
seem implausible that BHs accrete at the Eddington rate in 
halos in which the gas can be easily blown out after being 
photoionized. In this case, the Fitchett and Kidder models for 
the GW kick velocities can reach a maximum mass of only 
~ 2 X 1O''M0 by z = 6.43. In the upper right panel, we im- 
pose an additional requirement: halo mergers will lead to the 
mergers of their central BHs only if the mass ratio of the two 
halos, q, is between 0.1 and 1. Otherwise the merging halo 
is assumed to remain an orbiting satellite (with too long a dy- 
namical friction timescale for spiraling in), so its BH never 



reaches the center. In this case the Fitchett and Kidder mod- 
els reach a mass of only 3 x IQ^Mq at z = 6.43, and with the 
upper limit to GW kick velocities of Favata et al. the required 
mass of 4 X IQ'^Mq is barely reached, but for the Favata et al. 
lower limit to the GW kick velocities the required mass can 
still be reached at z ~ 9. 

We also show the result of raising ami„ to 50 kms~' in the 
lower right panel, and maintaining the requirement ^ > 0. 1 on 
the halo mass ratio. This higher limit might be more realistic 
if central BHs are only fed in galaxies with gravitationally 
unstable disks and a sufficient star formation rate (so that the 
gas does not all settle in a large, quiescent disk), and where the 
gravitational potential is sufficiently deep. Even in this case, 
the Favata et al. lower limit to the kick velocity can produce a 
BH mass above 4 x IO^Mq by z = 6.43, but with higher kick 
velocities the BH mass that can be reached is too low. Note 
that in the models of Haiman (2004J) BHs do not appear at 
all when a < a^in, whereas in our model BHs can still form 
and merge (but not accrete gas) for a < (Ji„j„. This, together 
with the higher v^sc we adopt, explains the difference between 
Haiman's results and ours. 

4. CONCLUSIONS 

Combining a continued mass growth of BHs from 
Eddington-limited gas accretion at high radiative efficiency 
(e ~ 0.1) with BH mergers, we can still account for the pres- 
ence of BHs with ~ 1O^M0 at z = 6.43, starting from stellar 
black holes of IOMq. Typically, these BHs can form by the 
merger of ~ 10^ stellar BH seeds, growing by a factor ^10'* 
in mass by gas accretion. We have shown that the ejection of 
some BHs due to GW recoil does not necessarily impede the 
growth of BHs in halo centers. 

Our model for the formation of these BHs involves two 
optimistic assumptions: First, that the BH masses grow by 
accretion with an e-folding time shorter than fso/ '^ 4 x 10^ 
years, for most of the time from their birth as stellar black 
holes at z '^ 20 to z — 6. Second, that BH mergers in which 
the recoil velocity is less than the halo escape velocity result 
in the BH immediately returning to the center and continu- 
ing its gas accretion, without significant losses from the en- 



ergy of the emitted GWs. These two assumptions are highly 
questionable. The first one may seem unlikely in view of the 
small fraction of galaxies that host quasars at lower redshift, 
although this fraction might be high in the earliest massive ha- 
los to form. Regarding the second, a more realistic calculation 
might show that black holes spend a lot of time finding their 
way back to the halo center after the GW kick, even when 
they do not escape the halo (" Hut & Reell 19921 iMerritt et all 
2004; Madau & Ouataert 2004). 

The earliest redshift at which a mass of ^ IO'^Mq can be 
reached strongly depends on the value of e assumed. Making 
e slightly greater than 0.1 would increase fso/ and reduce the 
number of e-folding times available between the epoch of for- 
mation of the first stellar black holes at z ~ 20 and the earliest 
time when luminous quasars are observed. This would then 
make it impossible to attain the required mass. At the same 
time, that is perfectly possible if black holes gain their mass in 
a short period of super-Eddington accretion, which is in fact 
the main alternative to the model of continuous Eddington- 
rate growth we have assumed here. 

In the absence of super-Eddington accretion, the scenario 
of continuous, luminous gas accretion and large numbers of 
mergers implies that quasars should be strongly correlated at 
high redshift. The halo hosting the luminous quasar at z = 6.43 
should have merged from lower mass halos in the recent past. 
For example, in the randomly generated merging history used 
for the realization shown in Figure 1, the halo mass at z = 9 is 
only about one third of its mass at z = 6.43, and several halos 
of mass greater than 10 "Mq merge with the main halo in the 
intervening time. Most of these halos would need to host their 
own quasar shining near the Eddington luminosity for most 
of the time, within a distance of the turn-around radius of the 
halo at z = 6.43 (a few comoving megaparsecs). 
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